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ji ! ABSTRACT 

6 : 

background environments can affect the star and planet formation processes oc- 
curring within them. In turn, young stellar members can shape the background 
environment and thereby provide a feedback mechanism. This work explores one 
\ aspect of stellar feedback by quantifying the background X-ray radiation fields 

produced by young stellar objects. Specifically, the distributions of X-ray lu- 
\ minosities and X-ray fluxes produced by cluster environments are constructed 

as a function of cluster membership size N. Composite flux distributions, for 

^ \ given distributions of cluster sizes N, are also constructed. The resulting dis- 

tributions are wide and the X-ray radiation fields are moderately intense, with 
the expected flux levels exceeding the cosmic and galactic X-ray backgrounds by 

X ; factors of ~ 10 - 1000 (for energies 0.2 - 15 keV). For circumstellar disks that 

c3 | are geometrically thin and optically thick, the X-ray flux from the background 

cluster dominates that provided by a typical central star in the outer disk where 
r > 9 — 14 AU. In addition, the expectation value of the ionization rate provided 
by the cluster X-ray background is (x ~ 8 x 10~ 17 s _1 , about 4-8 times larger 
than the canonical value of the ionization rate from cosmic rays. These elevated 
flux levels in clusters indicate that X-rays can affect ionization, chemistry, and 
heating in circumstellar disks and in the material between young stellar objects. 

Subject headings: stars: formation — planetary systems: formation — open 
clusters and associations: general 
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1. Introduction 

Most star formation in our galaxy occurs within embedded clusters, which are them- 
selves located inside giant molecular clouds. The radiation fields produced by these stellar 
nurseries can influence the formation of additional cluster members, and especially their ac- 
companying planetary systems (Adams 2010). More specifically, the processes of star and 
planet formation can be affected via [1] heating of starless cores, leading to evaporation and 
the loss of star forming potential (e.g., Gorti & Hollenbach 2002), [2] ionization within star- 
less cores, leading to greater coupling between the magnetic fields and gas (e.g., Shu 1992), 
thereby acting to suppress continued star formation, [3] evaporation of circumstellar disks, 
leading to loss of planet forming potential (e.g., Shu et al. 1993, Hollenbach et al. 1994, 
Storzer & Hollenbach 1999, Adams et al. 2004), and [4] ionization of circumstellar disks, 
which helps maintain the magneto- rotational instability (MRI), which in turn helps drive 
disk accretion (e.g., Balbus & Hawley 1991). With regard to circumstellar disks, we note 
that the background radiation from the cluster environment often dominates that produced 
by the central star, especially at ultraviolet (UV) wavelengths (e.g., Hollenbach et al. 2000; 
Armitage 2000; Adams & Myers 2001, Fatuzzo & Adams 2008); nonetheless, the photoevap- 
oration of disks is often driven by stellar radiation (e.g., Shu et al. 1993; Alexander et al. 
2004, 2005, 2006). 

Multi-wavelength observations of embedded clusters over the past two decades have 
provided a wealth of information on these environments (Lada & Lada 2003; Porras et al. 
2003; Allen et al. 2007), making it possible to assess what impact a given physical process can 
have on star and planet formation. Indeed, a comprehensive analysis of the effects of FUV 
and EUV background fields on cluster environments has already been performed (Armitage 
2000; Adams et al. 2006; Fatuzzo & Adams 2008; Hollenbach & Gorti 2009; Holden et al. 
2010). The main goal of this paper is to perform a similar analysis for the X-ray band. 

This paper is organized as follows. In Section 2, we discuss the cluster environment, 
present our characterization of the stellar IMF, and outline the basic approach for assigning 
X-ray luminosities to stars in our cluster models. We discuss the statistical aspects of our 
cluster model in Section 3 and present distributions of the X-ray luminosities for clusters with 
varying stellar membership size N. In Section 4, we construct the corresponding distributions 
of X-ray fluxes, both for varying N and for different forms for the stellar density profiles. 
The implications of these findings are discussed in Section 5; the paper concludes in Section 
6 with a summary of results and a discussion of potential applications. 
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2. The Cluster Environment 

Studies of clusters out to 2 kpc (Lada & Lada 2003) and out to 1 kpc (Porras et al. 
2003) indicate that in the solar neighborhood, the number of stars born in clusters with 
N members is (almost) evenly distributed logarithmically over the range N « 30 to 3000, 
with half of all stars belonging to clusters with N < 300. These stars are contained within 
a cluster radius R c ranging between 0.1 — 2 pc, with the radial sizes of observed clusters 
following an empirically determined law of the form 



R C (N) = lpc 



N 
300 



with 7 = 1/2 (see Figure 2 of Adams et al. 2006, which uses the data from Carpenter 2000 
and Lada & Lada 2003). For simplicity, we adopt this relation to specify the cluster radius 
throughout most of this paper, although we use 7 = 1/3 when considering a population of 
clusters that extend up to N = 10 5 in order to be consistent with the radial sizes observed 
for larger clusters (for further discussion and supporting data, see Chandar et al. 1999; 
Pfalzner 2009; Proszkow & Adams 2009). The gas within the cluster environment not used 
to build stars accounts for roughly 70 - 90% of the total (gas and stellar) mass of the cluster, 
and extends well beyond R c until it eventually merges smoothly into the background of the 
molecular cloud. 

We assume that the stars within a cluster are drawn from a parent population that is 
described by the universal initial mass function (IMF) observed in our Galaxy. The resulting 
probability distribution is therefore independent of membership size N. We adopt the broken 
power-law form of the IMF presented in Kroupa (2001; not corrected for binaries), expressed 
in terms of dimensionless mass m = M/M Q and truncated at m = 100, so that 

dN+ 



dm 

where 



Qm- a > } (2) 



a = 0.3 0.01<m<0.08, 

ai = 1.3 0.08 <m< 0.5, (3) 

a 2 = 2.3 0.5<m<100. 

The constants Cj are easily determined through the continuity of the IMF and the normal- 
ization of equation (j2J) to the number of stars within a cluster. For completeness, we note 
that the expectation value of the stellar mass for the adopted IMF is given by 

1 r 100 fdNA 
(m) = — m dm = 0.38 , (4) 



N+ Vo.oi V dm 
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where N* = J'(dN*/dm)dm. 

We now consider the X-ray luminosity of the stars within an embedded cluster. X-rays 
from low-mass stars (m < 2) are produced from flaring events in stellar coronae, whereas 
X-ray emission from high-mass stars (m > 18) is driven by stellar winds (intermediate mass 
stars most likely do not produce X-rays, with detections from such stars then attributed to 
low- mass companions). It is not surprising, therefore, that observations of X-ray luminosities 
of cluster members show a large variance in X-ray values (e.g, Feigelson et al. 1993, Feigelson 
et al. 2002, Preibisch et al. 2005); the variance also tends to decrease with the cluster age 
(see Alexander & Preibisch 2012). Nonetheless, these observations indicate that the X- 
ray luminosity depends on stellar mass. Modeling the exact nature of this dependence is 
hampered by the large variance in the data and by the fact that the derived correlation 
between Lx and m is model dependent (Preibisch 2005). In our work, we adopt a purely 
empirical approach and use a monotonic relation between stellar mass and X-ray luminosity. 
Specifically, we use the results obtained by Preibish et al. (2005), and focus primarily on 
the PS model for which a linear regression fit yields the relation 

log (L x [ergs -1 ]) = 30.34 + 1.13 logm , (5) 

and a standard deviation of A(logLx) = 0.64. As shown in Figure dj this best fit line (solid 
line) and the corresponding ±A(logLx) deviations (dotted lines) are in reasonable agreement 
with the low- mass data (squares) obtained from Chamaeleon I (Feigelson et al. 1993) and 
the low-mass data (gray area) and intermediated/high- mass data (circles) obtained from the 
Orion Nebula (Feigelson et al. 2002). For completeness, we also calculate cluster luminosity 
distributions using results from the SDF models, for which a linear regression fit yields the 
relation 

log (Lxfergs- 1 ]) = 30.37+ 1.44 logm , (6) 
and a standard deviation of A(logLx) = 0.65. 

As elaborated on in Section 3, we obtain X-ray luminosity distributions for clusters with 
N members by randomly sampling the IMF for each stellar member, and then assigning a 
luminosity to each star based on its mass. For a fixed stellar mass, we include a distribution of 
X-ray luminosities; guided by the aforementioned observations, this latter step is performed 
by randomly selecting a value of log [Lx{w)\ from a normal distribution with the mean given 
by equation (jSJ) or equation (jSJ) and with a standard deviation w = A(logLx). We note that 
for the PS model, the mean luminosity of the resulting distribution for stars of mass m is 

/oo 
P[x; 30.34+ 1.13 logm, 0.64]10 x cfe = 6.4 x 10 30 m L13 , (7) 
-oo 

where P[x;pL, w] is the probability density function of mean fi and standard deviation w. 
This value is denoted by the dashed line in Figure [TJ Since the values of log [L x {m)] are 
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selected from a normal distribution, the corresponding distribution of Lx values are not 
normal, but rather, are positively skewed. As a result, log Lx{m) > 30.34 + 1.13 log m (as 
seen in Figure 1). The expectation value of the X-ray luminosity — per star — over the full 
distribution of stellar masses is then given by 

1 Z" 100 - dN 
(L x ). = — L x {m)— -dm = 2.6 x 10 30 ergs" 1 . (8) 
N J om dm 



We next consider the contribution to the total X-ray luminosity that comes from stars 
with mass m for a parent population drawn from our adopted IMF. We start by noting that 
the function 

tpx{m) = L x (m) -j— , (9) 

which characterizes the X-ray luminosity produced by the stellar population (as a function 
of mass), peaks at m — 0.08 and drops off as a power-law ip x oc m -1 ' 17 for masses m > 0.5. 
As a result, most of the luminosity of a cluster per mass decade, which scales as mip x , is 
expected to come from stars with mass m > 0.5. To illustrate this point, we randomly 
select 10 7 stars from our IMF and then select the corresponding luminosity from equation 
(jHJ), including the variance about this mean relationship as prescribed above. The resulting 
X-ray luminosity distribution is shown in Figure [2J The distribution displays a broad peak 
near L x ~ 10 29 erg s -1 ~ Lx(0.08), and approaches a slope of —1 (as indicated by the 
dotted line) as Lx —> 10 31 erg s _1 ~ L~x(0.5). The distribution of X-ray luminosity per mass 
decade at higher values is therefore nearly flat, indicating that above m ~ 0.5, stars ranging 
in mass between m = 10 — > 100, although fewer in number, can contribute nearly the same 
to the total cluster X-ray luminosity as stars ranging in mass between m — 1 — > 10. This 
result is consistent with the weak dependence of the function mipx oc m~ a17 on mass above 
m = 0.5. It is worth noting that our calculated distribution, as shown in Figure 2, is in 
good agreement with the observed distribution of X-ray luminosities of the Orion Nebula 
population (see Figure 3a in Feigelson et al. 2005). 



3. X-Ray Luminosity Distributions 

With the IMF and the prescription for obtaining the luminosity of the sources specified, 
we now determine the characteristics of the X-ray luminosity distribution for a given cluster 
of size N, including both the expectation value and the variance. The X-ray luminosity for 
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Fig. 1. — Summary of observational data for the mass vs X-ray luminosity relation for stars 
in young cluster environments. Squares - low mass Chamaeleon I stars presented in Table 
6 of Feigelson et al. (1993); circles - intermediate/high mass ONC stars presented in Table 
5 of Feigelson et al. (2002); grey region - region of parameter space spanned by soft X-ray 
emission of low mass ONC stars, as presented in Figure 14 of Feigelson et al. (2002). The 
solid line shows the results of the linear regression fit obtained by Preibish et al. (2005) 
using the PS model, and the dotted lines mark the corresponding ±w deviations, where 
w = A(logLx) is the residual standard deviation. The dashed line represents the mean 
value Lx(m) as given by equation (JTj) in the text. 
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Fig. 2. — Probability distribution of the X-ray luminosities. This distribution was ob- 
tained by first randomly selecting a stellar mass from the IMF defined by equations (|2J) 
and 01]), and then obtaining the corresponding luminosity by randomly selecting Lx from 
a normal distribution of \ogLx with mean given by equation (jSJ) and standard deviation 
w = A(logLx) = 0.64 as discussed in the text. The dotted line has slope -1 and is included 
for illustrative purpose. 
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a single cluster, denoted as Lx(N), is given by the sum 

N 

L X {N) = Y, L *H (10) 

where L X j is the X-ray luminosity from the jth member, and N is the membership size of 
the cluster. In this formulation, we assume that the X-ray luminosity for a given star is 
determined by a distribution whose properties depend only on the stellar mass, and that 
the stellar mass is drawn independently from the specified stellar IMF. The sum in equation 
fTTU]) is thus the sum of random variables, where the variables are drawn from known and 
independent distributions (which are determined by the IMF and the Lx — rn relation). In 
the limit of large N, the expectation value (Lx)c of the X-ray power for the cluster (which 
is a function of N) is given by the expression 

(Lx)c = N(L x )*, (11) 

where (Lx)* is the expectation value of the X-ray power per star, as given by equation 
©• As usual, the central limit theorem implies that the distribution of values Lx{N) 
obtained by sampling over many clusters must approach a Gaussian form in the limit iV — > oo 
(e.g., Richtmyer 1978), although convergence is often slow. One of the issues of interest 
here is the value of stellar membership N required for these statistical considerations to 
be valid; similarly, we would like to know the fraction of the cluster population that has 
such sufficiently large N. In its limit of applicability, this Gaussian form for the composite 
distribution is independent of the form of the initial distributions, i.e., the shape of the 
distribution is independent of the stellar IMF and the mass-luminosity relation. The width 
of the distribution also converges to a known value given by the expression 

1 - 

( w )x = nY. w J (w) x = VNw 0} (12) 

3=1 

where w is the width of the individual distribution, i.e., 

wZ = (L x U-(L x }* 2 . (13) 

We can write wo = Q(L X )*, where Q is the dimensionless width of the X-ray luminosity 
distribution. Using the specifications given above, we find Q ~ 17. For a given membership 
size N, the mean value of the distribution of cluster luminosities converges to the value 
N(Lx}*, whereas the width of the distribution converges to VnQ(L x }*- As a result, the 
mean value will be larger than the width of the distribution only if the number of stellar 
members exceeds Nq = Q 2 ps 300. Small clusters (with N < Nq « 300) thus have such wide 
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distributions of X-ray luminosity that the mean value does not provide a good estimate. 
For larger clusters (with N > Nq), the X-ray luminosity is adequately characterized by the 
expectation value, but the distributions retain significant (relative) widths. These trends are 
illustrated in the figures that follow below. 

Figure [3] illustrates the X-ray luminosity distributions obtained from our cluster model, 
scaled by the cluster expectation value N(Lx)* from equation (TTTj) . for clusters with mem- 
bership size N = 10 2 , 10 3 , 10 4 , and 10 5 (short-dashed curves, from left to right). The 
long-dashed curve represents a Gaussian profile with mean \x = 1 and standard deviation 
w = Q/vTo 4 ", whereas the solid curve represents a Gaussian profile with mean /i = 1 and 
standard deviation w = Q/^10 5 . It is clear from this result that as N increases, the sam- 
pling becomes more complete, and the distributions approach a Gaussian form. However, 
in practice, even with large values of cluster membership size N, the asymptotic form is 
not reached. We note that the clusters in our solar neighborhood, as observed by Lada & 
Lada (2003) and Porras et al. (2003), are in the regime of incomplete sampling, and the 
distributions of their X-ray luminosities are expected to have median values that are smaller 
than N(Lx)* while retaining significant tails. 

Figure H] shows the X-ray luminosity distributions (with physical units) for clusters with 
stellar membership size iV = 100, 300, 1000, 3000 and 10000, constructed using both the PS 
model (solid curves) and the SDF model (dashed curves) for the X-ray luminosity vs mass 
relation. As expected, the mean/median values of the distributions increase with increasing 
N; in addition, the width of the distribution decreases with increasing N. The SDF models 
produce distributions with larger means and standard deviations (by about a factor of 2 
for larger clusters) than their PS model counterparts. These characteristics result from the 
stronger mass dependence for the luminosity in the SDF models, resulting in a broader range 
of luminosity values sampled in that case. 

4. X-Ray Flux Distributions 

In this section we construct distributions of the X-ray flux impinging upon the cluster 
members. Unlike the FUV and EUV bands, for which a significant fraction of the luminosity 
is produced by the most massive member of the cluster (e.g., Fatuzzo & Adams 2008), the 
total X-ray luminosity can arise from comparable contributions of numerous stars. As a 
result, the overall X-ray emission can be spread out over the entire cluster region, although 
this assumption is not universally true. For example, about 50% of the total X-ray luminosity 
from all (N > 1500) of the X-ray emitting stars in the Orion Nebula Cluster comes from 
the most massive star. Nevertheless, we can obtain a benchmark estimate for the X-ray flux 
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Fig. 3. — Distribution of X-ray luminosities. The short-dashed curves show the probability 
distributions for clusters with (from left to right) N = 10 2 stars, N = 10 3 stars, N = 10 4 stars, 
and N = 10 5 stars. Note that the luminosities (along the horizontal axis) are normalized 
by the average value one would obtain with complete sampling of the stellar IMF. For 
comparison, the long-dashed curve shows a Gaussian distribution with mean /i — 1 and 
standard deviation w = Q/y/10 4 , and the solid curve shows a Gaussian distribution with 
mean [i — 1 and standard deviation w — Q/ VlO 5 . 
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Fig. 4. — Probability distribution for X-ray luminosities. The five pairs of curves show the 
probability distributions for clusters of size (from left to right) N = 100, 300, 1000, 3000 and 
10000. Solid curves were obtained using the PS model; dashed curves were obtained using 
the SDF model. 
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within the cluster environment by dividing the expectation value of the luminosity N(Lx)* 
by the fiducial surface area A = nR c (N) 2 , where R C (N) is given by equation (TjQ), which yields 
a value of 

F ™ = ^7lk = 2 - 6 x 10 " 5 er § cm " 2 s ~* • ( 14 ) 

nR c (Ny 

Of course, the X-ray fluxes within a cluster, and within a collection of clusters, will 
take on a distribution of values, where Fxo represents an order-of-magnitude estimate for 
the mean. In order to construct these distributions of the X-ray flux, we first randomly 
select the position of a test star within a radius R C (N) (assuming a spherical distribution), 
weighted by a chosen density profile. For this study, we consider density profiles with the 
simple forms p oc 1/r and p oc 1/r 2 ; these profiles roughly bracket the possible profiles 
expected from observations. The probability that the test star is located a distance between 
r and r + dr from the cluster center is then given by 



pr 2 dr 



P(r) dr = — . (15) 



We then randomly selected the positions of the remaining N — 1 stars within the cluster, 
following the same prescription used to set the location of the test star. For each of these 
field stars, a luminosity is obtained as detailed in Section 2. The flux impinging on the test 
star is calculated, and the entire process is then repeated 10 7 times in order to build up 
flux distributions. Note that this procedure builds a slightly different distribution than the 
one that describes the flux values found within the cluster environment, for which case the 
location of the test star (which was weighted in accordance to the specified density profile) 
would have been replaced by a randomly selected location within the cluster that was not 
weighted by the density profile. 

The results are shown in Figure (for p ~ 1/r) and in Figure O (for p ~ 1/r 2 ). The 
range of fluxes realized within the clusters is comparable to (but generally less than) the 
benchmark estimate of equation ffl^l) . Note that the flux distributions peak at nearly the 
same locations for both density profiles (we use the same expression to set the cluster radius 
for both cases, so that the mean stellar density for a given cluster size N is the same for both 
density profiles). However, stars with a p oc 1/r 2 profile are more centrally localized, leading 
to a distribution that extends out to higher values of flux. As the stellar membership size 
N increases, sampling of the stellar IMF becomes more complete, and the flux distributions 
shift to higher values (they move to the right in the figures) and become somewhat narrower. 
In rough terms, however, the expected X-ray flux levels are given by Fx ~ 10~ 5 erg cm -2 
s -1 , with a factor of ~ 3 — 4 variation for p oc 1/r profiles and ~ 5 — 10 variations (above 
the peak) for p oc 1/r 2 profiles. 
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Fig. 5. — Probability distribution for the X-ray flux impinging upon stellar members for 
clusters where the stars follow a radial density distribution p oc 1/r. The five curves show 
the probability distributions (left to right) for clusters with varying membership sizes N = 
100, 300, 1000, 3000 and 10000. The vertical line denotes our benchmark value of X-ray flux 
defined by equation ( 1T31) . 
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We complete our analysis by constructing composite flux distributions for all stars 
within: i) a cluster population similar to that observed in our solar neighborhood, for which 
the largest cluster has N max = 3000 stars; and ii) a cluster population with cluster member- 
ship size extending up to N max = 10 5 . For case (i) we assume a cluster radius as specified 
by equation (0Q) with 7 = 1/2, whereas for case (ii) we assume 7 = 1/3. Following the 
observational result that clusters in the solar neighborhood have sizes N that are (almost) 
evenly distributed logarithmically, we assume that both populations are described by the 
cumulative distribution 

log(iV) - log(30) 
71 ' hg(N max ) - log(30) ' 1 ) 

where f(N) denotes the fraction of stars born in clusters of stellar membership less than or 
equal to iV (see, e.g., Figure 1 from Fatuzzo & Adams 2008). Note that for N max = 3000, 
/(300) = 0.5, so that half of the stars in the entire cluster population belong to clusters with 
N < 300; for comparison, this "median point" is much larger for the sample with N max = 10 5 , 
where /(1732) = 0.5. Composite distributions are then obtained by first sampling over the 
assumed cluster population to set the size of the cluster that our test star belongs to, and 
then performing the process outlined above to calculate the flux impinging on that star as 
a result of the other stellar members. This process is then repeated 10 7 times in order to 
build up a statistically valid distribution. The results are presented in Figures [7] and El 



5. Implications and Discussion 

With the X-ray flux distributions specified, we can now assess some of the corresponding 
effects on the cluster and on forming stars (and planetary systems). We first compare the X- 
ray flux received by circumstellar disks from the background cluster with that received from 
their central stars. If we ignore geometry, the sphere of influence of the central star (in X- 
rays) is defined by the radius where the X-ray flux provided by the background environment 
Fxc of the cluster is equal to that provided by the star Fx*, i.e., 

£ = °-° 39 -s ™ 2 s " (ws^ft) (srsur = F ™' (17) 

where we have scaled the X-ray luminosity to a typical value. For moderate-sized clusters 
with N < N max = 3000 such as those found in the solar neighborhood, and for p oc 1/r 
density profiles, the X-ray flux provided by the background cluster has a characteristic value 
Fxc ~2x 10 -5 erg cm -2 s _1 , but varies by factors of 3 — 4 on either side of this benchmark 
(see Figures [5] and Cj) . For clusters with p oc 1/r 2 density profiles, the X-ray flux has nearly 
the same peak value, but varies by factors of 5 — 10 (see Figure EJ). As a result, the sphere 
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Fig. 7. — Composite distribution for the X-ray flux impinging upon stellar members in a 
population of clusters similar to that observed in the solar neighborhood with N max = 3000. 
The solid curve was obtained by assuming a density profile of the form p oc 1/r. The dashed 
curve was obtained by assuming a density profile of the form p oc 1/r 2 . The vertical line 
denotes our benchmark value of flux as defined by equation ( |T4|) . 
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Fig. 8. — Same as Figure 7, but for a population of clusters with membership size extending 
up to N max = 10 5 , and with radii given by equation (CP where 7 = 1/3. 
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of influence of the star (in X-rays) extends out to about r « 1300 AU, much larger than the 
typical disk size. If we extend the cluster sample to larger membership sizes with N max = 10 5 , 
the characteristic value of the X-ray flux increases to Fxc ~ 6 x 10 -5 erg cm -2 s _1 (see Figure 
EJ) and the sphere of influence of the star has radius r « 760 AU, still somewhat larger than 
most circumstellar disks. 

Taken at face value, this finding would indicate that the X-ray flux for the planet-forming 
region of the disk is (usually) dominated by the central star rather than the background 
environment of the cluster. However, this comparison ignores geometrical factors arising 
because the disk is not perpendicular to the outgoing stellar radiation. The expected vertical 
extent Az of the disk, in the context of X-ray absorption, is estimated to be (Az)/r « 
0.3 — 0.4 (see Figure 1 of Gorti & Hollenbach 2009). As a result, the disks are relatively 
thin geometrically, and the intercepted X-ray flux will be significantly smaller than the 
value presented in equation (TIT)) . For the limiting case of a geometrically flat disk, the flux 
intercepted by the disk has the following radial dependence 

Mr) = £5^ «(!-«) bsi^bJ ' (18) 

where u = R*/r and where we have taken the large radius limit to obtain the final expression 
(Adams & Shu 1986). Using this expression, we find that the central star dominates the 
X-ray flux only out to r « 14 AU (9 AU) for clusters distributions with N max = 3000 (10 5 ). 
For these systems, the outer disk thus receives more X-ray flux from the cluster than from 
the central star. As an example, consider radial position r = 30 AU and the extended cluster 
distribution with N max = 10 5 ; the X-ray flux from the star (from equation [TB] ) is estimated 
as Fx ~2x 10~ 6 erg cm -2 s -1 , 30 times smaller than the background radiation flux. 

Although the background cluster can provide more X-ray flux to circumstellar disks than 
their central stars, the X-ray irradiation of the inner portions of the disks by coronal emission 
from young stars remains important. Several studies, complementary to the analysis of this 
paper, have addressed the process of photoevaporation driven by X-rays from the central 
star (e.g., Ercolano et al. 2008, 2009). The predicted mass loss rates can be substantial 
(typically of order 10~ 8 M Q yr _1 ) and most of the mass loss takes place at disk radii in the 
range 10 - 20 AU, near the transition where X-ray flux is dominated by the star (inside this 
radius) and the background cluster (outside). 

On the scale of the cluster itself, ionization of the cluster gas is an important considera- 
tion. The main contribution to ionization in star forming regions is often taken to be cosmic 
rays, which have a fiducial ionization rate ( C r ~ 10~ 17 s _1 (Shu 1992); more recent studies 
suggest a slightly higher value (cr ~ 2.6 x 10~ 17 s _1 (van der Tak & van Dishoeck 2000). 
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For comparison, the ionization rate due to the X-ray background radiation is given by 

( x = ^a(E x )^J(r,x ), (19) 
E x Ae 

where E x ~ 1 keV is the characteristic energy of the X-rays, a(E x ) is the total photoelectic 
cross section (see the discussion of Glassgold et al. 2000), the function J is an attenuation 
factor, which is, in turn, a function of the optical depth r in X-rays, and the threshold value 
xo = E Q /kT x (given by the lowest energy E of the X-rays under consideration and by the 
temperature T x of the emission region). The function J is expected to be of order unity 
(Krolik & Kallman 1983) and we will take J = 1 for simplicity. The remaining factor in 
equation (fl9|) is the factor E x /(Ae), where Ae ~ 35 eV is the energy required to make an 
ion pair. This ratio thus represents the enhancement factor resulting from individual X-rays 
making multiple ionizations (note that E x /(Ae) ~ 30). A power-law fit to the cross section 
results in the form 

a(E x ) = a (E x /lkeV)- n } (20) 

where the index n = 2.485 and the fiducial value of the cross section <jq = 2.27 x 10~ 22 cm 2 
(Glassgold et al. 1997). Putting all of the above results together, we find that to leading 
order the X-ray ionization rate becomes 

( x « 8 x 10- 17 s- 1 ( „* XC _ 2 zi ) • (21) 
\2 x 10 erg cm A s 1 / 

For the given background X-ray flux value, the X-ray ionization rate is larger than the 
cosmic ray ionization rate by a factor of ~ 4 — 8 (see also Lorenzani & Palla 2001). Keep 
in mind, however, that the cosmic ray ionization rate varies with galactocentric radius and 
with column density of the region (van der Tak & van Dishoeck 2000). Likewise, the X-ray 
fluxes have a broad distribution of values (as shown in Figures [5] - EJ) , so that the X-ray 
ionization rate can be much larger, or smaller, than the benchmark value given in equation 
( l2~Tj) . The relative importance of cosmic ray and X-rays for ionization is thus expected to 
vary significantly from cluster to cluster. 

The cluster is essentially optically thin to X-rays, so that young stellar objects within 
the cluster receive the unattenuated flux. We can estimate the optical depth as follows: The 
cross section of equation (|20|) can be integrated over the energy range of the observational 
surveys (the energy range used to construct the flux distributions) to find an average value 
(a) w (To/10 ~ 2 x 10~ 23 cm 2 . The number density n of a cluster with N = 100 - 1000 
and radius R c ~ 1 pc is about n = 3000 cm -3 . The mean free path of an X-ray is thus 
£ = I /nip) f« 5pc > R c . 

Another way to gauge the importance of the X-ray radiation fields is to compare the flux 
levels to those of the background galaxy and the background universe. Here we are interested 
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in X-rays in the range 0.2 - 15 keV, corresponding to the range in the observational surveys 
used to construct the X-ray luminosity vs mass relations used in this paper. The total flux in 
soft X-rays (in the range 0.28 - 1 keV) is estimated to lie in the range fx ~ 5 — 10 x 10~ 8 erg 
cm -2 s _1 (Fried et al. 1980), where we have assumed an isotropic radiation field to convert 
(energy-integrated) specific intensity to flux. This value has considerable uncertainty. In 
addition, comparable fluxes are present in harder X-rays, e.g., the range 1-15 keV that 
provides the other half of the range under consideration. These harder X-rays are thought to 
come from extragalactic sources (and are observed to be nearly isotropic), whereas the softer 
X-rays have a significant Galactic contribution (e.g., Boldt 1987). In any case, the total 
X-ray background flux is about fx ~ 1 — 2 x 10~ 7 erg cm -2 s _1 . Given the distributions of 
cluster properties considered here, and the distributions of fluxes for a given type of cluster, 
the X-ray flux provided by clusters F X c is larger than the background fx by factors in the 
range 10 - 1000. 



6. Conclusion 

This paper constructs the expected X-ray radiation fields provided by young embedded 
clusters. Specifically, we calculate the distributions of X-ray luminosity for clusters of a given 
membership size N and the corresponding distributions of X-ray flux. We also determine 
the distributions of X-ray flux for two cluster ensembles (for the cluster distribution found 
in the solar neighborhood and for an extended distribution with maximum membership size 
Nmax = 10 5 ). The flux distributions depend on the density profile of the stellar objects (the 
individual X-ray sources) within the cluster and we consider models that span the expected 
range of possibilities. The main results can be summarized as follows: 

The expected flux levels for X-rays are modest, with a characteristic value of order 
Fx ~ 1 — 6 x 10~ 5 erg cm" 2 s" 1 ; the range corresponds to the types of clusters under 
consideration, where the X-ray flux increases with cluster membership size N and with the 
central concentration of the density profile (Figures [5] -|S]). The distributions of X-ray flux 
are relatively wide, however, with the full width at half-maximum for log 10 Fx corresponding 
to a factor of ~ 3 — 4 for p oc 1/r density profiles and a factor of ~ 5 — 10 for p oc 1/r 2 
density profiles (see Figures |5] and [6]) . The cluster to cluster variation in X-ray flux will 
thus be significant. For completeness, we also note that X-ray emission is highly variable 
with time, and that the scatter in X-ray luminosity generally decreases with cluster age 
(Alexander & Preibisch 2012). 

Given the expected X-ray flux distributions for embedded clusters, as calculated herein, 
the outer regions of circumstellar disks (where planets form) receive comparable amounts 
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of flux from their central stars and from the background. The central star dominates for 
small clusters if one ignores geometric effects (equation |17]). whereas the background X- 
ray flux from the cluster dominates for geometrically thin disks (equation JT8] ) and also for 
larger clusters (Figure ED- Circumstellar disks are expected to be relatively thin and the flux 
impinging on them from the central star is significantly reduced by geometric effects. As a 
result, the background X-ray flux is expected to exceed that of the star for disk radii beyond 
r m 9 — 14 AU. We stress that significant variability in these values are expected given both 
the large variability in the X-ray luminosities of individual stars and the broad X-ray flux 
distributions within the cluster environment. 

Whereas disks receive comparable X-ray fluxes from their central stars and the back- 
ground cluster, the situation is markedly different for the case of both FUV and EUV radia- 
tion. At UV wavelengths, the background radiation fields of the cluster overwhelm those of 
the central stars for clusters with N > 100 (Armitage 2000; Adams & Myers 2001; Fatuzzo 
& Adams 2008). As a result, photoevaporation of circumstellar disks due to UV radiation 
is often dominated by the cluster background (not the central star), whereas mass loss due 
to X-rays can be dominated by the central star (and not the cluster). 

The ionization rate provided by X-rays in the cluster environment has a characteristic 
value Qx ~ 8 x 10 -17 s _1 (see equation [21]). This ionization rate is about 4-8 times the 
fiducial ionization rate (cr from cosmic rays, which are often considered to be the dominant 
source of ionization for molecular clouds and star formation (e.g., Shu 1992). The ionization 
rates for both X-rays (x and cosmic rays (cr have wide distributions, so that the relative 
importance of the two sources will vary significantly from cluster to cluster (and can vary 
with time). In a similar vein, the background of X-ray radiation produced by the cluster 
is larger than the background radiation provided by the galaxy or the universe. Over the 
energy range under consideration, the cluster contribution dominates by a factor in the range 
10 - 1000, with a typical value of ~ 100. 

The X-ray background radiation in young embedded clusters has a number of potential 
effects on star and planet formation (Glassgold et al. 2000; Feigelson et al. 2007). The ele- 
vated levels of ionization lead to greater coupling between the cluster gas and magnetic fields 
(e.g., Shu 1992); this increased coupling, in turn, leads to longer time scales for ambipolar 
diffusion and acts to slow down additional star formation. Enhanced ionization also affects 
MHD-driven turbulence, both in the cluster gas and in circumstellar disks via MRI (Balbus 
& Hawley 1991). The background X-ray flux will produce line diagnostics in circumstellar 
disks (e.g., Tsujimoto et al. 2005; Hollenbach & Gorti 2009) and will affect chemical re- 
actions in the gas (Aikawa & Herbst 1999). These processes, and many others, should be 
studied in the future to provide a more complete picture of how cluster environments shape 



-22- 



the solar systems forming within them. 
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